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ELEMENT  LOCATION  FOR  THE  ADA  ARRAY 

V.  C.  Anderson 

The  structure  of  the  hydrophone  element  mounts  used  in  the  ADA  array 
may  lead  to  slow  changes  in  the  actual  element  locations  during  the  operating 
period  and,  in  particular,  during  the  ascent  and  descent  operations  when  signi¬ 
ficant  changes  in  temperature  and  in  structural  stress  in  the  pressure  hull  can 
occur.  In  view  of  this,  we  intend  to  implement  an  acoustic  element  location 
system  which  can  monitor  the  positions  of  the  elements  in  nearly  real- time. 

The  approach  to  the  localization  system  is  to  establish  a  base-line  array 
of  high  frequency  reference  transducers  around  the  periphery  of  the  array.  This 
set  will  be  rigidly  coupled  to  the  deck  so  that  the  locations  can  be  measured 
during  installation  and  can  be  assumed  to  be  stable  within  the  required  tolerance. 
This  will  clearly  be  true  with  respect  to  distances  along  the  deck,  but  some 
curvature  or  deformation  could  occur  over  the  full  length  of  the  deck  in  what 
will  be  called  the  "z"  axis.  A  series  of  bootstrap  calibrations  down  the  full 
length  of  the  array  will  hopefully  measure  any  such  deformation.  If  the  cumu¬ 
lative  error  in  such  a  measurement  proves  to  be  too  large,  data  can  be  derived 
from  observation  of  a  distant  source  with  a  number  of  selected  dements  whose 
local  relative  geometry  has  been  determined  by  the  nearfield  calibration  system. 

The  high  frequency  base-line  array  will  be  used  to  survey  in  the  instan¬ 
taneous  locations  of  dual -frequency  calibration  sources  suspended  over  the  array 
on  the  inner  surface  of  the  protective  fabric  dome.  The  construction  of  these 
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sources  is  sketched  in  Fig.  1.  The  1  inch  cylinders  are  resonant  at  40  kHz, 
in  the  circumferential  mode  of  oscillation  and  couple  in  tiieir  interior  to  a  1/4 
lambda  resonant  chamber  operating  at  2.5  kHz.  They  are  pulsed  at  40  kHz 
for  localization  with  respect  to  the  base-line  array  and  then  are  operated  in 
the  continuous  mode  at  2.5  kHz.  For  phase  localization  of  the  elements  the 
sources  are  sequenced  in  time  so  that  phase  arrivals  from  the  different 
calibration  sources  can  be  measured  without  mutual  interference. 

The  following  analysis  indicates  that  a  set  of  6  transducers  in  the  base¬ 
line  array  and  6  dual  frequency  calibration  sources  will  suffice  to  measure  the 
element  locations  to  a  precision  of  +  -  .2  radians  at  2. 5  kHz.  This  will  limit 
tlie  degredation  of  the  array  pattern  on  axis  due  to  tlie  uncertainties  of  the 
element  locations  to  less  than  .  1  dB.  This  assumes  the  dominant  error  to 
be  in  the  phase  measurement  at  2.5  kHz  and  that  the  phase  measurement  error 
will  be  less  then  ±  .  1  radians. 

THE  LOCATION  ALGORITHM 

We  assume  three  fixed  points/ whose  positions  are  known,  the  velocity 
of  sound  to  be  known,  and  one  unknown  point  not  located  in  die  plane  of  die 
three  reference  points.  The  constants  in  the  problem  then  are  the  nine  coor¬ 
dinates  of  the  known  points  and  the  velocity  of  sound.  The  data  are  the  three 
travel  dmes  between  the  known  and  unknown  points.  The  basic  travel  time 
equations  are: 

2  2  2  2  2 
C  Tj  =  (x  -  Xj)  +  (y  -  y^)  +  (z  -  z^)  . 


(1) 
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using  a  Taylor's  expansion. 


(2) 


CT,2  »  (Xq  *  Xj)2  +  2(xQ  x.)  Ax  +  • . .  -KyQ  -  Yj)2  +  2(yQ  -  Yj)  Ay  +  •  •  •  +  <zQ  -  Zj)2  +  2(zQ  -  Zj)Az 

where  xQ  is  the  estimated  position  x  =  xQ  +  Ak  is  the  true  position. 

The  set  of  equations  can  be  expressed  in  matrix  form  as 

A  •  X  =  Y 


where 


X 

Vxr  yo ’ yr  Vzi 

(C2Tt  -  R*)/2 

y 

•  £  = 

xo“x2*  yo  " y2»  Vz2 

•  I  = 

(C2T2  -  R22)/2 

z 

X0-X3*y0-y3'Z0'Z3 

(C2T3  -  R32)/2 

and 


R)2=(xo‘xj)2  +  <y o*yJ,2  +  <W2 


solving  the  equation  for  X 


(3) 


X  =  A 


-1 


Y.  . 


(4) 


Cramer's  rule  is  used  for  the  matrix  inversion  rather  than  the  somewhat  more 
efficient  Gauss's  expansion  because  it  is  more  robust  in  an  automated  computation 
in  that  it  avoids  the  possibility  of  dividing  by  a  very  small  or  zero  coefficient, 
an  event  that  plays  havoc  with  the  dynamic  range  of  a  computer.  Cramer's 


rule  as  used  is: 


s ' |A|  k?i 


(5) 
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|  A  |  is  the  determinant  of  A,  A^  is  the  cofactor  of  <2^,  the  element  of  matrix 
A,  and  is  an  element  of  vector  Y.  The  determinant  Is  computed  using 
Laplace's  development, 

|a|  =#uAn  +  ^2A22  +  ^3A33  * 

In  the  computer  algorithm  the  magnitude  of  the  determinant  is  constrained  to 
be  larger  than  some  arbitrary  limit. 

Appendix  I  is  a  listing  of  the  subroutine  for  the  solution  using  Equation 
5  and  6.  A  modified  routine  is  listed  in  Appendix  II  which  constrains  the  maxi¬ 
mum  excursion  of  the  correcting  terms.  In  both  of  these  algorithms  a  generally 
slightly  higher  rate  of  convergence  is  achieved  by  the  use  of  a  gain  factor 
applied  to  die  correcting  terms. 

Convergences  for  several  cases  are  illustrated  in  Fig.  2.  The  upper 
curve  represents  the  extreme  case  where  the  estimate  is  placed  in  the  plane 
of  the  reference  transducers.  The  determinant  in  this  case  is,  of  course,  zero 
and  the  initial  excursion  is  finite  only  because  of  the  constraint  on  the  minimum 
size  of  the  determinant.  The  effect  of  the  additional  excurs!on  constraints  of 
II  is  apparent  in  the  central  cluster  of  curves.  The  number  of  iteration  steps 
is  reduced  by  30%  by  virtue  of  these  constraints.  However,  because  of  the 
additional  computation  required  for  the  constraint,  the  overall  processing  time 
is  actually  longer  for  program  II. 

The  gain  factor  of  1. 1  shows  a  significant  advantage  in  the  convergence 


of  these  central  curves,  but  not  for  the  lower  left  hand  ones  with  a  different 
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and  smaller  displacement  for  the  initial  estimate.  Here  the  gain  factor  actually 
reduces  the  convergence  rate  slightly. 

LOCATION  ERRORS 

The  statistics  of  the  errors  engendered  by  noise  or  instrumental  degra¬ 
dation  of  the  time  delay  data  can  be  determined  by  differentiating  Equation  5 


dV  “FJ/«T.dTl  + 


«F,/  _  dT.  +<*F./  ^  dT, 
1  aT2  2  j  oTj  3 


Assuming  the  dT's  to  be  independent  zero-mean  random  numbers  of  equal 
variance, 

crj2  =  ((  aFj/  orTj)2  +  (  aF  ^  aT^  +  (a^/  aT^  )  cr,2 

<rj/  crT  =  (1/  | A |  )((Aljj  Rj)2  +  <A*2J  Rj)2  +  (A^  r/ 

This  ratio,  the  ratio  of  the  location  error  in  the  x^  coordinate  to  the  time  delay 
error,  is  dependent  on  the  geometry  of  the  problem. 

Figures  3  and  4  display  error  ratio  contours  in  a  plane  of  height 
ZQ  =  .3,  0.3  times  the  reference  array  spacings.  From  these  contours  it  is 
apparent  that,  over  the  x-y  aperture  of  the  reference  array,  the  error  ratio 
in  the  xy  plane  at  a  height  of  .3,  corresponding  approximately  to  the  ADA  geometry. 
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w  11  be  less  than  2.0,  and  less  than  1.0  over  most  of  the  aperture.  The  z  error 
ratio  will  be  less  than  1.0  over  the  reference  array  aperture.  Notice  that  the 
error  ratio  increases  rapidly  for  locations  beyond  the  limits  of  the  aperture, 
the  z  error  ratio  increasing  much  more  rapidly  than  that  in  the  xy  plane.  The 
reason  is,  of  course,  that  z  is  being  measured  more  nearly  endfire  to  the  array. 
As  a  matter  of  fact,  it  is  the  coupling  of  the  z  error  into  the  xy  coord  nate 
estimates  that  gives  rise  to  the  rapid  increase  in  the  xy  error  beyond  the  ref¬ 
erence  aperture  boundary. 

For  comparison  purposes  the  error  contours  for  z^  =  1.0  are  shown 
in  Figs.  5  and  6. 

PHASE  LOCATION 

Once  the  dual -frequency  calibration  sources  have  been  located  using 
high  frequency  pulsed  measurement  of  travel  time,  they  are  operated  in  the 
low  frequency  continuous  mode  for  location  of  the  main  array  elements. 

In  this  case,  where  the  location  is  accomplished  by  the  use  of  steady  state 
sinusoidal  signals,  measuring  the  phase  at  an  array  element  from  a  succession 
of  source  locations,  die  above  algorithm  can  be  used  to  obtain  corrected  positions. 
One  way  of  accomplishing  this  is  to  convert  the  phase  measurement  to  a  time 
delay  error  as  follows. 

The  phase  from  the  j1*1  source  is  represented  by 

EXP  (i  </>  j)  =  COS  ( <^)  +  i  SIN  (  =  EXP  (i  <o  T.) 
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The  difference  between  the  phase  received  with  time  delay  and  that  which  would 
have  been  observed  for  the  range  to  the  estimated  position  is, 

BXP(I»T  -  =COS(<^)COS(^Rj)+SWj)SIN(^Rj)+i(SIN(^)COS(^R.)  -  COS(<f>.)SINUR 

The  difference  A  =  is  then,  within  an  ambiguity  of  ±  if  , 

A=  ft»T.  =SIN"1(SIN(«*»j)COSURj)  -  COS(<^)SINUR.)) 

and  also,  to  remove  ambiguities, 

A  =  o»T.  -  ifRj  =C08‘1  (COS(^)COSURi)+  SIN(^)SINUR^)) 

The  elements  of  the  vector  Y  of  Equation  4  can  be  approximated  more  directly, 
ifcTj-R^ZL  cTj  +  R^  by, 


2  2 
C  TJ.--RQJ 


“  VcVV'CR0i 


A 

(I) 


Subroutines  for  the  phase  rotation  and  conversion  of  phase  components 
into  time  delay  values  are  given  in  Appendix  3.  The  lookup  tables  used  in 
these  routines  are  sized  to  maintain  an  rms  error  of  less  than  .003  wavelengths. 
The  rotation  subroutine  also  includes  provision  for  the  correction  of  constant 
phase  shift  offsets  of  an  instrumental  nature. 

VELOCITY  OF  SOUND 

It  will  be  necessary  to  monitor  the  velocity  of  sound  in  the  water  inside 
of  the  dome  because  variations  with  temperature  can  be  significant  over  the 


MPL-U-43/76 
TM  279 


range  of  temperatures  which  may  be  experienced  In  operation.  Pairs  of  the 
base-line  reference  array  elements  may  be  used  for  this  purpose,  however, 
careful  attention  must  be  paid  to  their  location  to  assure  that  a  reliable  path 
exists  between  the  elements  in  the  horizontal  direction.  This  means  that, 
first  of  all,  the  elements  should  be  positioned  at  a  height  of  18  Inches  off  the 
deck  so  that  they  will  be  in  the  first  maximum  of  the  Lloyd  mirror  interference 
pattern.  Second,  there  should  be  no  main  array  elements  in  the  line-of -sight 
between  the  elements  to  be  used  for  velocity  of  sound  measurement. 

As  a  very  useful  check  on  overall  Instrumental  time  delay  offset  it 
will  be  helpful  to  install  an  additional  high  frequency  transducer  immediately 
adjacent  to  one  of  the  reference  elements. 

A  suggested  configuration  for  the  high  frequency  reference  array  trans¬ 
ducers  is  sketched  in  Fig.  7. 
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PEAL  CU jUUHD , OMEGA. 

:-m******MM****imS  ROUTINE  CONUERTS  CO  AND  QURD  PHASE 
%  COMPONENTS  TO  A  TIME  DELAV  MODULO  ONE  PERIOD'.  THE 

MAIN  PROGRAM  MUST  SET  UP  THE  LOOKUP  TABLE  WITH  STATEMENTS 
y.  such  ns . . . 

..  REAL  JiRRAV  ARCSINE C- 30: 503 ;  REAL  1 5 

y.  FCF;  i  :*  -50  STEP  1  UNTIL  50  DO  ARCSINECI3  »*  ARCS  IN  <  I /’SO.  >5 


BEGIN  V  •  ■ 

REAL  ‘NOOV TS  ‘  •  . 

MUD:=  SQRT<C0*#2+QUADif:s»i2>: 

CO:=CO*50/MOD;  QUAD: =QUAD*50/MOD* 

.707  THEN 
0  THEN 

ARCSINE CQUAD3  * 

ELSE  T  :=  SIGN<GUAD;>*<3. 1HI59-ABS<ARCSINEIQUAD3>> 
ELSE  7  :®  SI GN< QUAD1) #<1 .570 SO- ARCSINE IC01)S 
PHASETIME  :=  IOMEGA* 

END; 


IF  AbS<CO>  GEO 
IF  CO  GTR 
T  :  = 


100  PROCEDURE  ROTATE <CO»  QUAD*  RJ» KAY» SHIFT)  i  t  . 

101  REAL  CO» QUAD »RJ» KAY » SHIFT? 

102  ^♦♦♦♦♦♦♦♦♦♦ROTATION  OF  COMPONENTS  CO  AND  QUAD  BY  THE  PHASE  SHIFT 

103  %  OVER  PATH  RJ  AT  WAVENUMBER  KAY  AND  BY  INSTRUMENTAL  PHASE  SHIFT 

104  *  "SHIFT"  IN  RADIANS.  THE  MAIN  PROGRAM  MUST  SET  UP  A  SINE  LOOK 

1 05  *  UP  TABLE  AS  FOLLOWS -  • 

1 06  V.  REAL  ARRAY  SINE  CO:  1603  1  REAL  I*. 

107  y.  FOR  I  0  STEP  1  UNTIL  160  DO 

103  Y.  SINE CID  :*  SIN <3. 14159/1/64)  ) 

1 09  %♦♦♦♦♦♦♦♦♦♦ 

110  BEGIN 

111  REAL  ALPH . BET » CO 1 » QUAD 1 » 

112  INTEGER  TRUNC5 

113  TPUNC  :*  I NT  < <RJ*KflY+SHIFT> *20.3718)  » 

114  TPUNC  :*  TRUNC  -  TRUNC  MOD  1281 

115  ALPH  :=  SINEtTRUNC+323? 

116  BET  :«  SINECTRUNCDi 

117  CO 1  :■  CO*ALPH  +  QUAD*BET) 

118  QUAD1  :■  QUAD* ALPH  -  CO*BETI 

119  CO  :«  COll  QUAD  *•  QUAD1) 

120  END) 

ti 


DISTRIBUTION  LIST 


Director 

Advanced  Research  Projects  Agency 
1400  Wilson  Boulevard 

Arlington,  Virginia  22209  (3) 


Office  of  Naval  Research  Branch  Office 

1030  E.  Green  Street  ; 

Pasadena,  California  91101  (1) 


Mr.  T.  F.  Locas 
ONR  Resident  Representative 
University  of  California,  San  Diego 
P.O.  Box  109 

La  Jolla,  California  92037  (1) 


Cdr.  Donald  D.  Pizinger,  Code,  212 
Office  of  Naval  Research 
Department  of  the  Navy 

Arlington,  Virginia  22217  ...  .-.  (1) 


t  «  y ; * i ••  •  >.• : 


